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Abstract 

The DREAM collaboration has introduced the "Q/S Method" for obtaining 
the energy estimator from simultaneous Cherenkov and scintillator readouts of 
individual hadronic events. We show that the algorithm is equivalent to an 
elementary method. 
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1. Introduction 

The response of a hadronic calorimeter to an incident pion or jet of energy 
E can be written as 

hadronic response = E[f em + (1 - f em )(h/e)} (1) 

where a fraction f em is deposited in electromagnetic (EM) cascades, mostly 
initiated by 7r° decay gamma rays, and h/e is the energy-independent ratio 
of detection efficiencies for the hadronic and EM energy deposits^ (Here and 
elsewhere the energy E is normalized to the electron response.) In the case of a 
dual readout calorimeter, in which a Cherenkov signal Q and scintillator signal 
S are read out for each event, Eq. [I] can be generalized: [TJ [21 01 [5]. 

Q = E[f em + (l-f em )(h/e\ Q )] (2) 
S = E[f em + (1 - f em )(h/e\s)} (3) 

John Hauptman has suggested the less cumbersome notation rjx = (h/e\x), 
which we use in this paper: 

Q = E[f em + (l-f em )r) Q ] (4) 
S = E[f em + (1- f em ) vs ] (5) 



Email address: degroomSlbl.gov (Donald E. Groom) 
1 Whether one writes this ratio as e/h, as is conventional, or h/e is usually unimportant. 
This is not the case in Sec.pfl where we regard h/e as a stochastic variable. 
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The EM fraction f em is a feature of the event, while the efficiency ratios t]q and 
rjs are different for the two channels. Equations [4] and [5] are the starting point 
for any analysis of dual-readout hadron calorimetry data. 

If fem, VQj an d Vs are known exactly, and if there are no photoelectron 
or other statistical contributions, then Q and S are uniquely determined by 
the incident hadron energy E. If, on the other hand, all of these quantities 
are subject to statistical fluctuations, then E as determined from the equations 
must be regarded as the estimator of the hadron (or jet) energy for a particular 
event. 

These equations appear explicitly in Fig. II of the first DREAM paper at 
the Perugia Conference on Calorimetry in High Energy Physics [1] and appear 
either explicitly or implicitly in subsequent DREAM papers. The most complete 
description of the DREAM analysis is given by Akchurin, et al.[2] (henceforth 
Ak05), and it is the basic reference for this paper. Several data reduction 
schemes are described, but the algorithm considered most basic is the fairly 
convoluted "Energy-independent Q/S correction method." We show here that 
it can be obtained in a few lines from Eqns. [4]and[5j 



2. The energy-independent Q/S correction method. 

The estimator E, whose determination is the object of the analysis, can be 
eliminated by dividing Eq. [4] by Eq. [5j to obtain 

Q _ fem + (1 ~ fem)VQ 
S fem + (1 - fem)VS 

This is Eq. 2 in Ak05, except that in that paper values of h/e special to the 
DREAM experiment are inserted for tjq and rjs- It can be solved for f em . 
Although the result, 

, _ (Q/S)vs " VQ , 7 x 
Jem (l-r, Q )-(Q/S)(l- Vs ) ' [i) 

is not given in the paper, its availability is assumed in the rest of its discussion. 

Leakage corrections are incorporated as part of the Method. They are obvi- 
ously important, but here we assume they have already been made to Q and S 
as given m Eqns. g] and [5) 

The hnal estimator of the energy, called Sanai) is given by Ak05's Eq. 7: 



S'final — S c . 



(8) 



.1 + femPl/PO 

where pi/po = e/h — 1. We identify Sfi na i with the energy estimator E, and 
replace S COII by S because the leakage corrections hve already been made. From 
context, e/h is e/h\g- The equation then becomes 

F = o e / h \s (q] 

° 1 + f em (e/h\ s - 1) { > 

(10) 



VS + fem(l - T)s) 
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which we recognize as just a rearrangement of Eq. [5j 

It remains to insert the expression for f em into this equation. Simplification 
of the result is fairly tedious, but finally yields 



S 



(i - t,q) - (Q/s) V3 

Vs - VQ 



(11) 



3. Direct solution 

We can write the simultaneous equations [4] and [5] as 



S 

with immediate solutions 
E 

fem 



-(l-r) Q )\ (l/E\ 

-(I-VS)J [fern 



VQ 

Vs 



S 



(Hgh (Q/S)y s 

Vs ~ Vs 

(Q/s)y s - vq 
(i-vs)-(Q/S)(i- VQ ) 



(12) 

(13) 
(14) 



for the estimators of E and f em on an event-by-event basis. This method has 
been published elsewhere [U [5], but the identity of the approach with the "Q/S 
method" was not previously recognized. 



4. Discussion 



In part because of the relatively small number of particles involved early in a 
hadronic cascade, the efficiency with which the hadronic energy deposit is visible 
in either the Cherenkov or scintillator channel varies from event to event. In 
contrast, the efficiency with which the EM deposit is detected varies little. The 
result is that tjq and 77s are stochastic variables, mostly reflecting the variation 
of h. The values of tjq and 775 required to compute the energy estimator for each 
event via Eq. 11 (or Eq. 13) are not only unknown but unknowable, give "only" 



dual readout. In actual data reduction, there is little choice but to replace them 
by their mean values: 



E = S 



(i-(vq))-(Q/S) (vs) 



(Vs) ~ (vs) 
It is also useful to rewrite Eqs. |4]and[5j 

(Q/E) - f em + (1 - f em ) (vq) 
(S/E) = f em + (1 - f em ) ( vs ) 



(15) 



(16) 
(17) 



Since (Q/E) and (S/E) are linear in f em , (Q/E) is a linear function of (S/E), 
describing a line segment from ((Q/E) , (S/E)) — ((vq) , (vs)) at the all-hadronic 
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Figure 1: Energy-independent event locus in the Q/E—S/E plane. With increased energy, 
resolution improves and the mean moves upward along the locus. 



extreme, f em = 0, to ((Q/E) , (S/E)) = (1, 1), at the all-EM extreme, f em = 1. 
This event locus is shown in Fig. [I] As the energy increases, the Monte Carlo 
event scatter shown in the figure moves upward and becomes more clustered as 
the resolution improves. 

The energy-independent event locus has slope 

R ^-^l . (18) 

This slope can be determined either by linear fits to monoenergetic (test beam) 
event distributions in the Q/E-S/E plane, or, perhaps more accurately, by 
separately finding (tjq) and (77s) via n/e measurements as a function of energy. 



It can be used to cast Eq. 13 into a more tractable form[U[S]: 



Since f em is not needed in data reduction, it is only of academic interest. 
Experimental distributions based on DREAM data are shown in Refs. [1] and 
[2] (Ak05). These are broadened by resolution effects, and so do not necessarily 
conform to < f em < 1. 



Acknowledgment s 

John Hauptman's critical comments and suggestions have been particularly 
helpful. This work was supported by the U.S. Department of Energy under 



4 



Contract No. DE-AC02-05CH11231. 
References 

[1] R. Wigmans, "First results of the DREAM project," Proc. 11th Inter. Conf. 
on on Calorimetry in Particle Physics, Perugia, Italy, 29 March-02 April, 
2004, cd. C. Cccchi, P. Ccnci, P. Lumbrano, & M. Pepe, World Scientific 
(2005) 241. 

[2] N. Akchurin, et al., Nucl. Instr. and Meth. A 537 (2005) 537. 

[3] N. Akchurin, et al., Nucl. Instr. and Meth. A 550 (2005) 185. 

[4] D.E. Groom, Nucl. Instr. and Meth. A 572 (2007) 633. 

[5] J. Beringer et al. (Particle Data Group), Phys. Rev. D86 (2012) 010001. 



5 



